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We consider the excitation of large amplitude gyrotropic vortex core precession in a Permalloy 
nanodisk by the torques originating from the in-plane microwave current flowing along the inter¬ 
face of the Permalloy/BiaSea heterostructures, in which the huge charge-to-spin conversion ratio 
is observed |T]. We consider analytically and by micromagnetic modelling the dependence of this 
excitation on the frequency and magnitude of the microwave current. The analogies of the vortex 
dynamics and the Landau phase transitions theory is demonstrated. These findings open the pos¬ 
sibility to excite gyrotropic vortex motion with the current densities far lower than by any other 
means. 


Magnetic vortices are very interesting nanoscale mag¬ 
netic states which attracted research interest from both 
fundamental and applied perspectives J2HI5]. Magnetic 
vortices are formed by in-plane magnetization that curls 
clockwise or counter-clockwise around the small volume 
of out of plane magnetization in the center which is 
called the vortex core. This magnetization distribution in 
nanopillars mm represents an energy favourable state 
for a wide range of geometries [HI [19]. 

The increased attention to magnetic vortices is largely 
related to their excited states, which can be maintained 
by an external stimuli (field, current etc.). Among 
the dynamic modes of magnetic vortices in nanopillars, 
the gyrotropic mode is of the greatest interest, because 
it is the lowest energy collective excitation B I unl¬ 
it is represented by the circular (or gyrotropic) mo¬ 
tion of the vortex core, accompanied by corresponding 
high-amplitude oscillations of the mean magnetization of 
the whole nanopillar. There are promising reports on 
the possibility to build vortex-based non-volatile mem¬ 
ories [22 Spin Transfer Nano-Oscillators (STNOs) 
mm, spin diode frequency sensors [22], and even more 
sophisticated systems based on collective vortex dynam¬ 
ics [24] . 

The key issue in this context is the manner of the vor¬ 
tex gyrotropic mode excitation. Initially, it was proposed 
to excite vortices by an external microwave magnetic 
field (see for instance Refs. [12 122 [25]). An alterna¬ 
tive approach is to use dc spin-polarized current to ex¬ 
cite the gyrotropic motion of the vortex (see for instance 
Refs, m m H])- Although the latter approach is eas¬ 
ier to be implemented in practice and appears to be more 
energy efficient, the issue of reducing the required current 
amplitude remains a high priority. 

An original solution could be to use pure spin currents, 
which can lead to the generation of large spin torques act¬ 
ing on magnetization. Such torques, which can be due 
to the interfacial Rashba effect but more often are the 


result of the bulk spin Hall effect, are considered for the 
ferromagnetic metal layer with asymmetric Pt and A10 x 
interfaces [27] as well as for a graphene/Pt interface [281 . 
Also such torques are observed in the case of ferromag- 
netic/Pt [23] or ferromagnetic/Ta [30] bilayers, and more 
generally in systems with a ferromagnet in contact with a 
layer of large spin orbit coupling material (Pt, Ir, W, Ta, 
Pd etc.). Using these torques, it is possible to excite do¬ 
main wall motion [22132], to prevent Walker breakdown 
133] , to switch the magnetization of a ferromagnetic disk 
[301 [541 :55i and also to excite magnetization oscillations 
[3511371. In this context, an important question is to find 
new materials with properties favourable the generation 
of larger spin currents. 

Nowadays, bismuth selenide and bismuth telluride ma¬ 
terials with giant spin-orbit interaction (SOI) attract 
a lot of attention mainly in the contest of topologi¬ 
cal insulators (TI). The possibilities of magnetization 
switching [35] and enhancement of the Walker break¬ 
down [39] by TI were demonstrated theoretically. In 
2014, the bias current induced spin polarization in Bi 2 Se 3 
was observed experimentally [40] . Moreover, the mag¬ 
netization switching through giant spin-orbit torque in¬ 
duced by an in-plane current in a chromium-doped TI 
bilayer heterostructure was demonstrated experimentally 
at low temperature |41j . In [42] spin-orbit torques in 
Bi 2 Se 3 /Co 4 oFe 4 oB 2 o devices were evaluated at different 
temperatures using ferromagnetic resonance measure¬ 
ments. It was demonstrated that the magnitudes of the 
spin-orbit torques decrease be factor of ten with increase 
of the temperature from 50 I\ to 300 K , which indicates 
that origin of the observed spin-orbit torques in this sys¬ 
tem are topological surface states (TSS). Most of these 
appealing TSS-based effects, however, were observed at 
low temperatures and for very pure materials. Consid¬ 
erably ’’dirty” bismuth selenides and bismuth tellurides, 
which can be or doped either have many structure de¬ 
fects, still have a giant SOI, and thus are very promising 
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for spintronics. In this case the Fermi level is located in 
the conduction band, therefore, because of a very small 
statistical weight of TSS, the spin-orbit effects are deter¬ 
mined mainly by the band states. Recent experimental 
works demonstrate that charge current flowing at the sur¬ 
face of a thin film of Bi 2 Se 3 at room temperature can ex¬ 
ert a spin-transfer torque on an adjacent ferromagnetic 
permalloy (Py) film, with the strength greater than of 
any other sources of spin-transfer torque measured thus 
far fl] @3]. Moreover, in case of doped material (like 
BiTel) it was demonstrated experientially that although 
it ceases to be an insulator and demonstrates bulk con¬ 
ductance, it still possess strong spin-orbit coupling [44] . 
which expected to be enough to induce the magnetization 
dynamics. 

This is one of our objectives: to evaluate the feasi¬ 
bility and effectiveness to excite gyrotropic vortex mo¬ 
tion using torque originating at the BUSes/Py interface. 
The studied system (see inset in Fig. [T]) is represented 
by a NigiFeig disk with radius Rd = 100 nm and vor¬ 
tex magnetization distribution on top of a Bi 2 Se 3 film. 
Such a system is very attractive for direct observations 
of magnetization dynamics by optical means since the 
ferromagnetic layer is open at the top in contrast to con¬ 
ventional STNOs or spin-diode effect which require a cap¬ 
ping electrode [36]. First, we demonstrate by micromag- 
netic modelling that although the torque associated with 
the direct in-plane current leads only to a finite vortex 
position shift, the ones associated with radio-frequency 
(rf) in-plane current along the interface can efficiently 
excite large rotation of the vortex core. We consider the 
dependence of this excitation on the frequency and the 
magnitude of the rf current. The analytical description 
of this phenomenon is proposed. 
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Figure 1. (Color online) Vortex core orbit R CO re as a function 
of frequency of the input RF signal for jo = j so (black line), 
jo = 1/2 x j BO (green line) and jo = 3 xj so (red line), where j so 
is the characteristic current density for the considered system. 
Inset: the system under investigation from afar. The black 
arrows represent magnetization direction, while the red region 
represents the vortex core with positive m z . 

Here we study the magnetization dynamics in the nan¬ 


odisk under the action of an rf in-plane current. The disk 
thickness is d = 8 nm. The magnetization dynamics in 
the NiFe disk is described by the Landau-Lifshitz-Gilbert 
(LLG) equation with an additional term responsible for 
the spin-transfer torque caused by TI: 

M = - 7 MxH e// + -^(MxM)+T S TT, (1) 
M s 

where M is the magnetization vector, 7 is the gyromag- 
netic ratio, a is the Gilbert damping constant, Mg is 
the saturation magnetization and H e ff is the effective 
field consisting of the magnetostatic field, the exchange 
field and the anisotropy field. The spin-transfer torque 
T stt can generally be represented by two components: 
a parallel (in-plane) torque T|| = 7 T||m x ( e x x m) and 
a perpendicular torque Tq = 7 Tj_(m x e^,), where e^, 
is a unit in-plane vector perpendicular to current direc¬ 
tion and m is a unit vector along the magnetization. In 
materials such as Bi 2 Se 3 , a giant Rashba interaction is 
commonly present. This interaction causes significant 
spin accumulation at the surface [ 10 ], that can interact 
directly by exchange with the magnetic moment of the 
adjacent NiFe disk. The diffusive spin current inside the 
adjacent ferronragnet (similar to the spin Hall effect) can 
create a parallel torque, and the combined actions of the 
Rashba effect and the exchange interaction can produce 
a perpendicular one. The microscopic details of these 
effect in studied system are still not completely clear, 
but the first experimental evidence of the existence of 
both torques has already been obtained pQ. In this case, 
the torque densities can be expressed as t so = — divjg, 
where jg is a spin current. Substituting the spin torque 
ratio 9 = (js/j){2e/h), one can obtain torque density 
t/° = h0ij/2ed (with i =[|,J_), where j is the current 
density and e > 0 is the charge of the electron. 

It needs no mention that the spin-transfer torque T stt 
in general case may have more complex form than the 
one considered here (see Ref. [15]). However for the 
magnetic vortex, which is formed by the in-plane curl¬ 
ing magnetization and a tiny out-of-plane core, the dy¬ 
namical properties are defined solely by the in-plane vol¬ 
ume. Consequently in our case we can consider only the 
torques action on the in-plane component of the magne¬ 
tization, without any loss of generality. In other words 
m(r) = (cosx(r), siny(r), 0 ) T , where m(r) is a unit vec¬ 
tor along the magnetization and y(r) is the angle be¬ 
tween X-axis and magnetization direction. Using it, gen¬ 
eral spin-transfer torque form from Ref. [45] is reduced to 
Tq = TitmxeJ and Ty = T/mx (e^. xm) respectively, 
when one consider the vortex-based system. 

To investigate the possibility of vortex core excitation 
by the current flowing along the interface Bi 2 Se 3 /Py, 
we have performed a series of simulations using our 
micromagnetic finite-difference code SpinPM based on 
the fourth-order Runge-Kutta method with an adaptive 
timestep control for the time integration and a mesh 
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size 2 x 2 x 10 nm 3 . The larger mesh size in the z- 
direction has been chosen because the modelling the mag¬ 
netization almost does not change with the z-axis. To 
focus on TI action, we have not included the Oersted 
field in our simulations, because its value is not enough 
for vortex core expulsion and even for significant exci¬ 
tation. The Py magnetic parameters used in the mod¬ 
elling are: Ms = 800 emu/cm 3 , the exchange constant is 
A = 1.3 x 1CU 6 erg/cm, a = 0.01 and the bulk anisotropy 
is neglected. 

We find in micromagnetic modelling that the in-plane 
direct current cannot excite vortex gyration, but only 
slightly shifts the vortex core from the equilibrium posi¬ 
tion. Thus, subsequently, we will focus on the influence of 
rf currents. According to the micromagnetic modelling, 
after the application of an in-plane rf current with a den¬ 
sity of j = jo sin(u;£), the vortex starts gyrotropic motion 
after a transitional period. Then it reaches a station¬ 
ary orbit due to the action of the spin torque created at 
the Tl/Py interface. The dependences of the vortex core 
orbit Rcore on the frequency of the input RF signal for 
three different alternating current density amplitudes jo, 
are shown in Fig. [l] with the vortex core being resonantly 
excited. At the same time, an increase of the current am¬ 
plitude jo leads to an increase in the stationary orbit of 
the vortex core. Moreover, if the value jo becomes higher 
than critical value jg r , the vortex core gets ejected out 
of the dot (corresponds to the region on Fig. [l] where the 
red line overtakes the Rd value). It is also worth men¬ 
tioning that the excitation has nonlinear behaviour for 
high current values, similar to the Ref. m ■ For further 
analysis, let us introduce the dimensionless currents ra¬ 
tio j/jso, where j so is the characteristic current density 
for the considered system, and it is chosen in order to 
introduce the figure in arbitrary units. 

For analytical insight into the vortex excitation we an¬ 
alyze the Thiele equation which can be deduced from 
Eq.Q: 

G(e z x R) = fc(R)R + DR-F ST (2) 

where R is the vortex core position, G = —2npMsh/'y 
is the gyroconstant, p is the core polarity, and e z is the 
unit vector along the z-axis. The confining force is given 
with fc(R) = ujoG(1 + aR 2 /R 2 D ) where the gyrotropic 
frequency is ojq = ^ 7 Msh/Rn and a « 0.25. The 
damping coefficient is D = aG{\ ln(^ a ) + |) [3B] where 

l e = J ~ mentioned before, the spin-transfer 

force Fst consist of two contributions: a perpendicu¬ 
lar torque and a parallel torque. The first one modi¬ 
fies the energy of the vortex with the additional term 
E _l = — <5(H e ff x e z )R, therefore it can be represented 
as Fj_ = 6(e z x H e ff), where H e ff is the effective field 
corresponding to the perpendicular torque. This formula 
is already in use in the case of external magnetic field 
[5]. Because the perpendicular torque has exactly the 


same symmetry, hence it can be used here with the effec¬ 
tive field: H e ff = t ]° sin(aA)e x . Micromagnetic mod¬ 
elling gives S = 5.6 . Using the Feldtkeller 

ansatz [47, for describing the vortex configuration, the 
in-plane torque contribution can be represented by F y = 

/3G(cos(0)e x —sin(</>)e y ), where /3 = p c V^ri/ 0 sin(wt), 

p c is the vortex core width, V is the vorticity and </> is 
the polar angle of the core. Using this, Eq. 0 can be 
represented in polar coordinates (R, <f>) in the following 
form: 


R 

R 

<t>-- 



k(R) 

~G~ 


S H{t) 
G~R~ 

DR 
~ GR 


cos(</>) + 4 sin( 2 </>) 

R 

5 H(t) . . . /3 , 

+ -q—r- sm ( 0 ) + ^ cos( 2 </>) 


(3) 

(4) 


where H(t) = rf° sin(<*;f). 

The dependence of the stationary orbit of the vortex 
core on the amplitude jo of the rf current density is rep¬ 
resented in Fig. [2] The resonant frequency of excitation 
increases with the increase of the current amplitude jo 
(see inset in Fig. [2]b). Hence, we consider two different 
dependencies of the orbit on the amplitude jo : with fixed 
frequency, which is resonant for the average value of jo 
(see Fig. §a), and with adjusting of the frequency for 
every jo (see Fig. [2]b). The first regime is more easily 
implemented, while the second one demonstrates greater 
efficiency of excitation, because it always deals with res¬ 
onant frequency. The dependencies of the orbit R core on 
the amplitude jo for constant frequency, which is reso¬ 
nant for jo = jso, obtained by the micromagnetic mod¬ 
elling (black dots) and by numerical solution of Eq. §■© 
(red line), are shown on Fig. [2ja. The jump of the func¬ 
tion, corresponding to a sharp change of orbit, arises from 
the fact that at the left side of the current j so , which is 
resonant for the selected frequency, the excitation effi¬ 
ciency is reduced due to the frequency mismatch. At 
the same time, this frequency mismatch is compensated 
by increasing of the amplitude at the right side of the 
resonance current. 

In Fig. [2]b, the dependence of the orbit on the ampli¬ 
tude jo for the frequency is adjusted for every jo obtained 
by the micromagnetic modelling (blue crosses) and by the 
numerical solution of Eq. 0.0 (red line). It should be 
noted that the orbit in case of fixed frequency is always 
less than in case of adjusted frequency, except for the 
jump point, where they are equal. As follows from these 
results, the analytical solution demonstrates good coin¬ 
cidence with the micromagnetic predictions. 

Both the numerical simulations and analytical Thiele 
modelling demonstrate that the instant frequency of vor¬ 
tex core rotation oscillates during the period with the de¬ 
viation amplitude less than 4%, while the mean frequency 
is equal to the frequency of the alternating current. Tak¬ 
ing this into account, we can integrate Eq. 0.© over 
the period under the approximation that the frequency 
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Figure 2. (Color online) Vortex core stationary orbit versus amplitude jo of the rf current density, a) The orbit obtained 
by the micromagnetic modelling (black dots) and by numerical solution of Eq. M (red line) in case of constant frequency, 
which is resonant for the average value of jo = jso- b) The orbit obtained by the micromagnetic modelling (blue crosses), by 
solutions of Thiele equations Eq. 00 (red line), and by Landau like theory solution of Eq. 0 (dash black line) in case, when 
the frequency is adjusted for every current amplitude jo- The quantity j so here is chosen in order to introduce the figure in 
arbitrary units. Inset: the dependence of the resonant frequency of excitation on the current amplitude jo- 


is constant. In this case, some oscillating terms vanish, 
while others return constant values after integration, and 
the Eq. 0 then takes the form: 

p = —Dloqp — Dojoap 3 + h (5) 

where p = R/Rd, D = D/G and h = /2GRd- Let 
us define function R as p = —dR/dp. Then using Eq. 0 
function R can be represented in the form: 

R = -hp + Duj 0 ^~ + Duj 0 a^- ( 6 ) 

This function is similar to a typical free energy functional 
in Landau theory, if we consider p as an order parameter 
and h as an external field conjugate to the order parame¬ 
ter. In this notation, Eq. 0 is equivalent to the Landau- 
Khalatnikov equation. The stationary orbit due to Eq. 0 
can be found by solving OR/dp = 0. This condition coin¬ 
cides with energy functional minimum condition. These 
analogies are valid not only for rf torques produced by 
TI, but also for all possible oscillating torques with the 
same symmetry (for example produced by the Rashba 
interaction or by the external magnetic field or by the 
Oersted field). In this case, the dependence of the sta¬ 
tionary orbit of the vortex core on the amplitude jo of 
the rf current can be represented in the form: 

h/Dcuo = p + ap 3 (7) 

The solution of Eq.([7]) is shown on Fig. [2jb (dash black 
line). It demonstrates good agreement with Thiele mod¬ 
elling and the numerical solution of Eq. 0 .@ until the 
orbits about Rd/2 = 50 nm. For orbits greater than this 
value, the larger oscillations of the frequency during the 


period appears, which makes our approximation no more 
valid. 

Experimental observation of spin torque ratio for 
Bi 2 Se 3 with an 8 nm permalloy layer atop p,j gives 
6 1 | > 3.5. It is worthy to note that this value already 
takes into account shunting by the Permalloy layer. To 
avoid shunting problems completely and therefore to ob¬ 
tain the maximum possible efficiency, the permalloy layer 
can be changed by the insulating ferromagnetic layer, as 
it proposed in Ref. [lj. At the same time, the latest re¬ 
sults of the same group [45] give lower estimations for 
the considered system, which are about ss 6 ± ss 1 . 
Using these values of spin-torque ratio let us consider 
the magnitudes of the torques with field-like symmetry, 
which is responsible for vortex gyration considered in 
this work, in different cases. In case of conventional 
perpendicular current injection (CPP) the magnitude 
T ± P = where £ = 0.1 - 0.4 and P is the cur¬ 

rent polarization for a Magnetic Tunnel Junction with 
typical values P = 0.3 —0.4. In case of excitation by spin- 
orbit effects the magnitude r [ 0 = Therefore, the 

comparison between CPP and spin-orbit excitation gives 
r i°/ r i PP = jp ~ 10 in case of Bi 2 Se 3 . These estimations 
give us the following value of j so = 1.5 x 10 6 A/cm 2 . 
It means that in Permalloy/Bi 2 Se 3 -like heterostructures 
the created torque can be significantly more effectively 
than in case of perpendicular current injection, which is 
mostly used nowadays. 

In conclusion, we have demonstrated the possibility of 
vortex oscillations excitation by the torques caused by 
Bi 2 Se 3 -like structures. Using micromagnetic modelling, 
we prove that in-plane rf current along the interface of 
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the Bi 2 Se 3 /Py heterostructure can excite vortex core ro¬ 
tation with a required current density several times lower 
than previously observed. We consider the dependence 
of vortex excitation on the frequency and the magnitude 
of the rf current both by micromagnetic modelling and 
by analytical theory. An analytical description of the 
system under consideration is proposed. The analogy 
between the used approach and the Landau theory of 
phase transitions is demonstrated. On this basis, the 
vortex excitation by the RF current through the the 
interface of the B^Sea/Py heterostructures becomes a 
very promising option for practical applications in spin- 
transfer nanooscillators, spin-torque diodes and so forth. 
This new fundamental mechanism of vortex excitation 
and peculiarities of its dynamics remain valid not only for 
considered material, but also for any alternative sources 
of the spin-currents with in-plane spin polarization, like 
Bi-based materials and other materials with strong spin- 
orbit coupling. However, other materials likely to provide 
smaller magnitude of the effect than Bi 2 Se 3 . 
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